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It is widely recognized that event-by-event analysis of experimental data on mul- 
tiparticle reactions (and especially studies of fluctuation patterns seen there) provides 
us with very important and sensitive tool in our attempts to understand dynamics of 
heavy ion collisions They are particularly useful in searching for some special fea- 
tures of the quark-gluon plasma (QGP) equation of state 0. The question of the best 
method of their investigation allowing for the most information to be gathered is there- 
fore of great interest (cf., for example, reviews and references therein). Some time 
ago, a novel method of investigation of fluctuations in even-by-event analysis of high 
energy multiparticle production data was proposed and applied to nuclear collisions 
0. It is based on a suitably defined measure $, 



N 



where Z = Zi, (1) 



=1 



which originally was supposed to discriminate wheather or not fluctuations of a given 
observable x is exactly the same for nucleon-nucleon and nucleus- nucleus collisions 
Here where x denotes the mean value of the observable x calculated for all 

particles from all events (the so called inclusive mean) and is the number of particles 
analysed in the event. In (|l]) {N) and (Z^) are averages of event-by-event observables 
over all events, whereas the last term is the square root of the second moment of the in- 
clusive z distribution. By construction $a; = for independently produced particles @]. 

However, application of this method is not free from controversy. When first applied 
to NA49 data for central Ph — Ph collisions at 158 A-GeV |0 it apparently revealed 
that fluctuations of transverse momentum (x = pt) decreased signiflcantly with re- 
spect to elementary NN collisions. This in turn has been interpreted as a possible 
sign of equilibration taking place in heavy ion collisions, providing thus an enviroment 
for the possible creation of QGP. It was immediately realised that existing models of 
multiparticle production are leading in that matter to conflicting statements [0]. The 
more recent NA49 data |^ reported, however, a new value (almost an order of magni- 
tude greater then the previous one), which was the one corresponding to a pion gas in 
global equilibrium A number of attempts followed, trying to clarify the meaning 
of $ (cf., for example, [0, 11 1 and references therein.) In the mean time, it was ex- 



tended to study event-by-event fluctuations of " chemical" (particle type) composition 
of produced secondaries |jl2[, to study azimuthal correlations among them (which are 



important for studies of flow patterns observed in heavy ion collisions |T^) and to 
cover also higher order correlations [l^. Finally, $ has been also analysed by means of 



the nonextensive statistic both for correlations [ffH] and for fluctuations of chemical 



composition as well 
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In |]T6[ we have pointed that, if there are some additional fluctuations (not arising 
from quantum statistics, hke those caused by the experimental errors), which add in 
the same way to both terms in deflnition (|l]) of it would perhaps be better to use 
another form of it, like for example 



{N) 



(2) 



where the variances cr^ would cancel (being present only implicitly). We would like to 
elaborate on this criticism towards $ (as well as $*) measure in more detail here. Our 
point is (cf. also 0) that inclusive experiments provide us both with single particle 
distributions P{x) (i.e., with information on fluctuations of the x-values) and with 
multiplicity distributions P{N) (i.e., with information on fluctuations of and, when 
put together with P{x), also with information on correlations between x- value and 
multiplicity A^). Also correlations between produced particles (i.e., between x- values), 
especially Bose-Einstein correlations (BEG) resulting from their statistics, are to a 



large extend known ^7j. Measures $ (or $*) depend on all of them, because (notice 



that = 0"? = al 
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{N)al + {N{N - 1)) ■ cov(x„Xj) + c(x, A^), 



(3) 



where covariance cov(xj,Xj) 



pa"^, i.e., is given in terms of the variance 
a'^ of the variable x and their mutual correlation coefficient p, while the last term 
describes the possible correlation between the variable x and the multiplicity A^. Notice 
that assuming, for example, (x) n for given A^ being given by 

N-{Ny 



{x) 



N 



X 



1 + a 



{N) 



one gets for correlation term 



c{x,N) = ■ x^ ■ [{N^) - (N)^) . 



(4) 



(5) 



Therefore $ = is achieved only for independent x-values, in which case they are also 
uncorrelated. If this is not the case they are given by: 



ai + 



{NjN-l)) 
(N) 



■ COV {Xi,Xj) + 



{m) - {NY 



■a'^x^ 



{N)-l + 



a 



N 



a 



(N), 



■ cov(xi, Xj) + 



$($ + 2a,). 



(6) 



(7) 



3 



Because in reality (neither in experiments nor in models or event generators at- 
tempted to describe them) conditions of independence are not fulfilled, $, $* 7^ and 
a question arises: what information do they convey? In particular, can we learn more 
from event-by-event analysis performing it by means of the $ (or $*) measure rather 
than by being satisfied with the known inclusive distributions alone? In this respect 
one should notice that 0: 

• Both $ and $* depend strongly on the multiplicity fiuctuations and this de- 
pendence vanishes only in the limiting case of a Poisson distribution (for which 
{N^} - {N}^ = {N) and (A^(A^ - 1)) = (A^)^) where, for example, 

$* = (AT) ■cov{xi,Xj) + o?x^. (8) 

• Experimental data allow us to study the dependence of (x) at on the multiplicity 
N in an explicit way. In the case when such correlations are absent, the 
measure reduces simply to a two-particle correlation measure 

= ((iV)-l + ^ ■cov(a;„a;,). (9) 

• The infiuence of the two-particle correlation term cov(a;j, Xj) depends on the kind 
of particles involved: $ = for Boltzman statistics, $ < for fermions and $ > 
for bosons. However, this type of correlations, especially for BEC for bosons, 
are subject to extensive experimental and theoretical investigations where one 
mesures or modells the so called 2— particle correlation function C2{xi,Xj) [0. 
Because 

cov{xi,Xj) = [ [ dxidxj Xi^— Xj-^— [C2{xi,Xj) — 1], (10) 

it means that $ and measures are nothing but the correlation measure C2 
averaged over single particle distributions dn/dx (i.e., we are in fact not gaining 
but losing some information contained in C2 and dn/dx and none of $'s provide 
us with any new information). 



The above remarks, being obvious and essentially known, 0, 10 , have so far not 



been backed by any convincing numerical illustration. By this we mean calculations us- 
ing an event generation algorithm which would satisfy all conservation laws (especially 
energy-momentum conservation) and at the same time model also the BEC. Because 
both points seem to play a major role in a proper description of $ (or $*) measure 
[|], ^, which was not checked properly |21|], we would like to fill this gap by calculat- 
ing $ and $* in simple models of hadronization using an algorithm which preserves 
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both energy-momentum conservation and the original single-particle distributions and 
models at the same time all features of BEC in such processes |jl8|. To this end we 
shall compare $ and $* measures as function of mean multiplicity {N^~^) of negatively 
charged particles for hadronization processes proceeding without and with BEC. Two 
simple models of hadronizations of mass M will be considered: the cascade model 
(CAS) developed by us recently [0 (where the whole space-time and phase-space 
history of the hadronization process is explicitly known) and simple statistical model 
(MaxEnt) based on information theory approach proposed in (where details of 
hadronization are not available; both models were already used in ||18||). 



Most of the above mentioned applications of $ measure concerned transverse mo- 
menta X = pt 0, ^, 0, which is always positive. The corresponding variable here 
would be \p\ results for which are shown in Fig. 1. We have found it interesting to 
enlarge analysis for the variable which can take any sign, as momentum p in our case, 
results for which are shown in Fig. 2. The differences are striking. Whereas behaviour 
of $(x = \p\) is similar to that for x = pt, that of $(x = p) is much more sensitive to 
the limits imposed by the energy-momentum conservation (both measures are negative 
here). In both Figs. 1 and 2 we present $ and $* for CAS and MaxEnt models with- 
out and with BEC, in this later case for two different choices of the weights specifying 
BEC [^: constant P = 0.5 and Gaussian P as discussed in (which lead to very 
different BEC patterns). The results are given as functions of the mean multiplicity 
{N^~^) of negatively charged secondaries produced in hadronization process of mass M. 



The special feature emerging from our calculation is the fact that the effect caused 
by BEC is clearly visible in all cases (it is maximal for the case of constant weights 
because the BEC is maximal there [0). To make dependence on BEC more clear we 
present in Fig. 3 for x = |p| case the corresponding differences 5$ = ^BEC—^noBEC and 
= — ^noBEC- "^^^ curves show the best power fits of the type 5$ oc {N^~^)^ 

exhibiting very strong and growing dependence of the effect of BEC on (N^^^). There 
are some other characteristic features of the results worth attention, namely: {i) the 
scale of effect is different in $ and $* this reflects differences in their definitions (^) 
and (^; (ii) their different dependencies on {N^~^) indicate that ax plays a significant 
role here; [in) the completely different behaviours for x = p and x = \p\ (compare Fig. 
1 with Fig. 2) illustrate how changes dramatically whether the measured variable 
is restricted to being positive or not. Because in the second case the role of the energy- 
momentum constraints is much more important, dependence on this constraint is more 
pronounced. Notice that points for MaxEnt in Figs. 1 and 2 are consistently higher 
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than for CAS. This should probably be attributed to different ways of imposing conser- 
vation laws in both models: they are satisfied globally in the statistical model MaxEnt 
and locally, in every branching point, in the cascade model CAS. Because of this CAS 
prefers symmetric distributions in momenta (here forward-backward) whereas MaxEnt 
allows more frequently for asymmetric ones (for example, when one particle with large 
momentum is balanced by a bunch of low momenta particles in other hemisphere). The 
difference between both types of models seen in Fig. 3 originate from different shapes 
of C2{pi,Pj) functions, which are broader (what corresponds to smaller "radius" R) for 
CAS type of models. 

To conclude: we have demonstrated using numerical algorithm which preserves 
energy-momentum conservation and at the same time models BEC in multiparticle 
hadronization that $ (and also $*) measure is very sensitive both to the constraints 
provided by conservation laws and to the effects of correlations (exemplified here by the 
BEC). Both features must therefore be carefully accounted for when attempts are made 
to reach some conclusions concerning new physical effects when using these measures 
to the event-by-event analysis of data. We have also written down exphcit relations 
between both measures and all other observables (well defined and known in statistics) 
available from inclusive experiments. They show exphcitly that both $'s are closely 
related to inclusive correlation functions. In this respect we believe that no new in- 
formation is obtained from these measures in comparison with what is available from 
inclusive measurements. 
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Figure Captions: 



Fig. 1 $ (in GeV, left panels) and (in GeV^, right panels) as function of mean 
(negatively charged) multiplicity {N^~^) calculated for a; = |p| for two different 
hadronization models, CAS and MaxEnt, without and with BEC. Upper panels 
contain BEC obtained using constant weights, lower panels contain BEC cor- 
responding to Gaussian weights |2^. Black symbols denote CAS and MaxEnt 



models with BEC, open symbols without BEC. Stars correspond to MaxEnt and 
circles to CAS models, respectively. Lines are just interpolating between calcu- 
lated points. 

Fig. 2 The same as in Fig. 2 but this time for x = p. 

Fig. 3 Differences between events with and without BEC for x = \p\, 5$ (upper part) 
and 5$* (lower part) for models and BEC weights presented in Fig. 1 (with 
the same meaning of different symbols). Circles denote CAS and stars MaxEnt 
models with, respectively, constant (for open symbols) and Gaussian (for full 
ones) weights for BEC. Lines indicate attempts of best power-like fits of type 
(5$ oc (N^^^)^ with S equal to (going from top) to 2.0, 2.0, 1.5, 0.7 for upper 
panel and 3.6, 3.6, 3.1, 2.2 for lower panel. 
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